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1. Introduction 

Recent experimental results indicate that the quark-gluon plasma produced in the 
heavy- ion collisions stays strongly coupled at temperatures above deconfinement [||]. 
Therefore understanding the nature of QCD matter at high temperature and density 
requires non-perturbative techniques. Lattice QCD, being an intrinsically Euclidean 
formulation is not well-suited for calculating certain important dynamical observables 
such as transport coefficients or any sort of real-time correlation functions. 

For this reason holographic techniques based on the AdS/CFT correspondence 
[0], have recently attracted much attention in study of such dynamical phenomena. 
For example the shear viscosity [Q], jet quenching parameter |^ and the drag force 
[^] has been calculated with better success than corresponding perturbative findings. 

One approach for constructing gravitational backgrounds dual to QCD-like the- 
ories, is to search for deformations of ten dimensional AdS^ x 5*^ background that 
breaks supersymmetry and conformality. Such models enjoyed success in repro- 
ducing certain IR phenomena but they also bear some non-realistic features such as 
presence of KK modes arising from the extra dimensions. Another, more phenomeno- 
logical approach 0, instead of attempting at deriving QCD from fundamentals of 
lOD critical string theory, aims at deriving a 5D gravitational background from the 
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basic requirements of QCD. This idea goes under the name of AdS/QCD^] and also 
achieved partial success, especially in the meson sector. 

Generally, finite temperature in the holographic approach is introduced by com- 
pactifying the Euclidean time direction with period 1/T. One such obvious solution 
is the thermal graviton gas. Other more non-trivial solutions involve black-holes. 
The black-hole solutions correspond to the deconfined phase of the corresponding 
gauge theory [^, hence encode physics above the deconfinement transition. The 
purpose of this paper is to review the thermodynamics in the aforementioned 5D 
models. 

In the next section, we review the thermodynamic properties of the AdS/QCD 
models based on the hard-wall (HW) and the soft-wall (SW) geometries. By ex- 
tending the analysis of |jlO|], we compute quantities such as the energy, entropy and 
speed of sound as functions of T and compare them with the expectations from the 
lattice. Unlike the HW, the SW model shows good agreement with the lattice data. 
However, being a non-dynamical model the black-holes in this geometry do not obey 
the laws of BH thermodynamics. From a practical point of view this fact renders 
computation of certain quantities like the bulk viscosity ill-defined. Moreover, it does 
not give insight in the nature of the deconfinement transition. 

In section 3, we study a dynamical model based on dilaton-gravity [JTT], [T^ which 
is close in many respects to real QCD. This model is based on general expectation 
of stringy holographic QCD whose thermodynamic properties were derived in [|l^ 



We show that this background solves most of the problems in the AdS/QCD models 
at once, sheds light on the role of the gluon condensate in the phase transition and 
yields very good agreement with the lattice data. In the final section, we summarize 
the results, and discuss further directions. 



2. Themodynamics of the AdS/QCD Models 
2.1 Hard- wall model 

The simplest 5D holographic model for QCD is introduced in P]. The idea is based 
on the fact that QCD behaves nearly scale invariant for a wide range of energies 
ranging from far UV down to medium energy scales. Thus the authors of proposed 
a geometrical set-up based on the 5D AdS space with a cut-off in the deep-interior 
of the holographic coordinate r. The cut-off is introduced in order to break the 
conformal invariance in the IR, and eventually to model color confinement. The 
location of the cut-off at r = tq, is dual to the dynamically generated energy scale of 
QCD as Aqcd ~ ^/ro- We shall refer to this solution, as the "AdS cavity" for short. 

The model captures many basic features of QCD: one finds a discrete glueball 
spectrum by studying the fluctuations of the metric, an area law for the Wilson loop 
by studying classical string embeddings [Q, etc. However, the real success of the 
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model is in the meson sector and indeed the intention of the authors of [Q was to 
apply it there. In all of the 5D models that are discussed in this paper, the meson 
sector is generated by space-filling D4 and D4 branes. The fluctuations of the brane 
fields produce the meson spectra. In addition to reproducing certain generic features 
such as chiral symmetry breaking, existence of Nambu-Goldstone fields, Gell-Mann- 
Oaks-Renner relation, one finds %9 agreement with experimental data in the 1^, 0^ 
and spectra. 

Yet, the model is crude in many ways, especially when it comes to the glue sector. 
Running of the gauge coupling is not taken into account; not only the electric but also 
the magnetic quarks are confined (the 't Hooft loop also exhibits an area law); there 
is an ambiguity in computation of the glueball masses and related to this there is a 
degeneracy in the 2~^^ and 0^"*" glueballs |ff2| ; both the glueball and meson spectra are 



quadratic for large orbital quantum numbers. This list can be largely expanded, but 
here we shall focus our attention on the thermodynamics of the model and show that 
the model does not correctly fulfill expectations for the finite temperature physics 
either. 

2.2 Thermodynamics of the Hard- Wall 



Some thermodynamical aspects of the hard- wall model is investigated in |I0|. At 
finite temperature there are two competing solutions with the same asymptotics on 
the boundary: i) a thermal graviton gas which is just the AdS cavity with compact 
Euclidean time of circumference 1/T ii) the AdS black- hole solution with horizon at 
r = Th- The temperature is related to the location of horizon as T = l/vrr/i. As 
one heats up the black-hole, the horizon expands and at a particular temperature it 
coincides with the IR cut-off = tq. This is the minimum temperature for presence 
of the black-hole inside the AdS cavity Tmm = l/^rro. To find out the true minimum 
of the free energy, one computes the action evaluated on i) and ii) and then one takes 
the difference, see iTDI for details. 



Let us define the IR scale as A = 1/ro. Then, the free energy density^ for the 
hard- wall model, for T > Tmm reads^. 



fHw = {M/fA' 



(2.1) 



Here Mp is the Planck scale and i is the AdS radius. One finds a confinement- 
deconfinement phase transition at Tc = 2^A/7i. The pressure density is given by 
Phw = —fHW- 



word on notation: Wc shall define the thermodynamic densities as the thermodynamic 
function divided by the volume of the 3D space V3 times the number of degrees of freedom N^. For 
example the entropy density is s = S/{V^N'^). 

^Our notation for the action reads S = —MpN^ J y/g{R + V). Note that this involves an extra 

factor of 2 with respect to llQl where his k is related to our Mp as k"^ = Mp N^. 
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In order to compare the analytic results of the hard-wall model with the lattice 
data for QCD, one should fix the various energy scales in the model, i.e. Mpi and 
A. The latter is usually fixed by comparing the vector meson spectrum of the model 
with the lattice data0. One obtains A = 323 MeV. This, in particular yields a 
transition temperature at Tc = 122.3 MeV [|T0|, (see eq. (pTT|)). 

Fixing the Planck mass is more tricky. The most rigorous way is by comparing 
the high T asymptotics of QCD and the holographic model. At very high tempera- 
tures, the (quenched) QCD becomes a free gas of gluons with a limit value for the 
pressure density Pqcd/T^ 7r^/45 as T — > oo. In the hard-wall model, we find from 
( |2.1|) that the same quantity limits to (Mp£)^7r^. Equating the two yields, 

M/=(457r2)-i (2.2) 

We stress that this is a universal, model independent way of fixing the Planck mass: 
One obtains the same value for all of the models discussed in this paper. This is 
guaranteed to happen quite generally, if the geometry asymptotes to an AdS black- 
hole near the boundary. 




Figure 1: Comparison of the energy, entropy and pressure densities in the HW model 
with the lattice data of Boyd et al. (dashed curves). 



Having fixed all of the parameters in the model, we can compare the thermo- 



dynamic functions derived from the HW model with the lattice results. From ( |2.1| ) 
follow all thermodynamic quantities by standard rules. The entropy density can ei- 
ther be found by s = —df/dT or by the Bekenstein-Hawking formula which relates 
it to the area of the horizon. Eventually, the energy density follows from e = / + sT. 
All in all, one has. 



SHw = 4(M/)V^T=^ 



3.m'^ 



(2.3) 



where Mpl is given by ( p.2|) . Energy, entropy and pressure are compared in fig. |l[ 
Clearly, there is poor agreement. In particular s/T^ is constant in the model as a 
result of the underlying AdS geometry and e/T^ is a decreasing function unlike in 
QCD. 
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The latent heat is defined as the energy density at the phase transition. The 
lattice value[|l^] is Lh = (0.77Tc)^. From (|2.3|) one finds a finite latent heat also in the 



HW model, = 8{Mpi)^A^ = (0.97Tc)^. As this is a finite quantity, the transition 
is of first order. 

Presence of a first order deconfinement transition is in accord with our expecta- 
tions from large Nc QCD |T^. However, there are other shortcomings of the model. 
First, the conformal anomaly is a non-trivial function of T in QCD, see fig. ^ 
whereas this functional dependence is lost in the hard-wall model. One can compute 
this as / {y2,Nl) = chw ~ '^Phw from ( p?TD and (|2.3| ) and one finds a constant 



= SN^V^^MpiY . Similarly, the speed of sound can be computed as = s/cy 
where = de/dT is the specific heat of the system, and one finds = 1/3. This, 
of course refiects the fact that the underlying geometry is AdS, and is in complete 
disagreement with QCD where Cg is again expected to be a non-trivial function of T, 
see fig. 0. 

Secondly, when one computes the bulk viscosity from the Kubo's formula (see 
[0 for a recent treatment) one finds that (/s = which is again in disagreement 
with QCD. This latter result is rather disappointing because (/s is considered to be 
an important observable probing the quark-gluon plasma at RHIC, and its profile 
as a function of T reveals important information regarding the nature of the phase 
transition. In particular, both from the low energy theorems and lattice studies 0], 
it is expected to make a peak near Tc. 

Although there are numerous shortcomings of the the HW model^, it should be 
viewed as a first step in a holographic approach, instead of a rigorous construction. 
Indeed, even the fact that such a simple model captures certain basic aspects of QCD 
is astonishing and should be taken as a starting point for a deeper investigation. 

2.3 Soft- wall model 

Motivated by the partial success of the HW modelin the meson sector, Karch et al. 
introduced an improvement in [^] that softens the breaking of conformality in the 



IR. This is achieved by replacing the hard-wall at Tq by a non-trivial dilaton profile, 

0(r) = (Ar)2. (2.4) 

This introduces a dimensionful parameter A that sets the scale of the problem in the 
IR. The geometry is still taken to be AdS^. 



^The consistency of the model is also questionable. As discussed in and motivated by the 
critical string-theory constructions such as fl^ , the IR brane in the AdS cavity is viewed as an 
"end of space-time" as opposed to a boundary. However, from the point of 5D Einstein gravity, 
the IR brane really acts as a boundary of the geometry in the deep interior. Thus, in principle 
one should allow for a Gibbons-Hawking term also at the location of the IR brane. The authors of 
p7| investigated this issue and found that the deconfinement transition goes away, once a Gibbons- 
Hawking term is added at the IR brane. 
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Figure 2: Comparison of the energy, entropy and pressure densities in the SW model with 
the lattice data of Boyd et al. The diamonds refer to lattice data. 



One nice feature of the soft- wall model is linear confinement: The meson spec- 
trum is linear for large orbital excitation number and for large spin, as opposed to the 
quadratic spectrum of the hard-wall|^. However, some of the unphysical features 
in the HW carry over in the glue sector. In particular, there is no running gauge 
coupling^, and magnetic quarks are confined. 




5 T, 




T 

5 Tc 



Figure 3: Comparison of the trace anomaly and the speed of light in the SW model with 
the lattice data of Boyd et al. (dashed curves). 



Another main issue is that, the model is non-dynamical, i.e. it does not fol- 
low from a 5D gravitational action. Instead, the metric and the dilaton profile are 
imposed by hand^. Related to this, computation of the glueball spectra from the 
bulk-fluctuations is ill defined. Below, we shall see other problematic features at 
finite T^. 

^One my think of <i>(r) as a holographic dual to a running coupling, but this identification is 
problematic for non-dynamical fields. 

^In a dynamical Einstein-dilaton-tachyon theory is constructed that admits SW as a solution. 
However it is hard to understand the presence of Tachyon both in the gauge theory and in gravity. 

^The authors of |Q did not intend to apply the model to the glue sector. As a phenomenological 
model designed to describe the meson physics in the quenched approximation it is indeed appropriate 
and the question of whether it solves the equations of motion is not crucial. As we discuss below, 
it becomes crucial when applied to thermodynamics of glue. 



6 



2.4 Thermodynamics of the Soft-Wall 



The study of thermodynamics on this background is initiated in Once again, 
one considers two competing solutions at finite T: (i) the SW geometry with compact 
Euclidean time, (ii) AdS black-hole, appended with the non-trivial dilaton profile 



( p.4|) . As already mentioned, the construction is non-dynamical, hence neither of 
these two geometries solve the equations of motion of a 5D Einstein-dilaton system. 
One assumes that they are solutions to some unspecified gravitational theory and 
computes the free energy density with the prescription described in section ^]2| [|1^ 
The result is, 



fsw = 2{M/fT' 



(2.5) 



Here A is the parameter that appears in (|2.4| ) and Ei is the exponential-integral 



function. One obtains a phase transition at = 0.4917A. As before, one can fix 
the value of A by matching the lowest p meson mass and one finds, A = 338 MeV 
which yields a better than HW, Tc = 191 MeV |Q. The latent heat also turns 
out better than the HW model. One finds Lh = (0.725Tc)'^ that is very close to the 
Lucini et al.'s lattice result of (0.77Tc). 

The non-dynamical feature of the model manifests itself in the computation 
of entropy. The entropy as computed from the Bekenstein-Hawking formula and 
from s = —df/dT above do not match. The BH geometry does not obey the laws 
of thermodynamics, which makes the findings questionable. However, let us press 
on, and assume that one indeed obtains a free energy of the form (|2.5| ) from some 
unspecified dynamical theory and work out other thermodynamic functions. The 
computation is just as in the HW case and the results are summarized in figs. and 
^. These results are in very good agreement with the lattice study of P3| . 

It is surprising that a non-dynamical theory, constructed with many assumptions 
yield such good results and it begs for a better understanding. We shall, in the next 
section, investigate a dynamical dilaton-Einstein system with solutions similar to the 
form ( p.4|) in the large r region. 

An underlying dynamical theory is needed also to compute certain important 
observables such as the bulk viscosity (. In the holographic set-up, this quantity is 



computed using Kubo's formula |T6|. The reason this computation is ill-defined in 
the SW model is that one needs to solve for the bulk fiuctuations in an holographic 
computation and this requires that the background solves Einstein's eqs. 

3. Non-critical holographic QCD 
3.1 Dynamical Models 

There is a long history of the dilaton-gravity systems in the context of the AdS / CFT 
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correspondence. Due to lack of space, we are not able to provide an exhaustive list 
of references here. Rather we shall mention a few articles that are closely related to 



our approach. The papers |^ demonstrated that type string theory provides a 
fruitful set-up for gravity duals of running gauge coupling. They considered a lOD 
background that involves a dilaton and a bulk tachyon field. Asymptotics of the 
dilaton in the deep interior exhibits a log-running of the gauge coupling! (however, 
presence of the tachyon is confusing as there are no obvious dual gauge invariant in 
the gauge theory). Similarly, Gubser analyzed a dilaton flow in the context of 
type IIB, truncated to 5D''. Other notable papers that study a dynamical dilaton 
flow in the 5D set-up are 
similar to ours^ 



, [^, [T^ and pSf. The latter uses an approach very 

recently analyzed the 



Finally, Gubser and collaborators |2^, [ 

dilaton-gravity system at flnite temperature, obtaining results that are quite similar 

to m. 



3.2 Improved holographic QCD 

There are many reasons supporting a non-fermionic (such as type string theories) 
and a non-critical holographic approach From an economic point of view, flve 



dimensions provide all the necessary degrees of freedom to construct a dual of QCD: 
four dimensions where the gauge theory lives, plus a radial direction dual to the 
energy scale of the gauge theory. Furthermore, a brief study of the low energy degrees 
of freedom of 5D non-critical string theory yields a nice correspondence between the 
various objects in string theory and gauge theory [0. Absence of extra dimensions, 
hence absence of the undesired Kaluza-Klein degrees of freedom is another attractive 
feature. 

The only non-trivial bulk fields required to model the low energy dynamics of 
large Nc QCD are the metric (dual the the energy-momentum tensor), the dilaton 
(dual to Xym and TrF^) and the axion^ (dual to 9ym and TrF A F). Here, we shall 
present such a set-up []TT], |T2| and describe its zero temperature solutions. A simple 
5D action is, 



^5 = -M'pN^ 



i?- -(9$)2 + y($) 



I d^xVh K. (3.1) 

JdM 



where ^ is a yet undetermined potential for the dilaton. The second term above is 
the Gibbons-Hawking term, K being the extrinsic curvature on the boundary^^. 



^An early work on dilato n fl ow in the IIB set-up is 
^See however below eq. (3^) for various differences. 

^We will not be concerned with the axion in this paper. Its action is suppressed in the large Nc 
limit, hence can be consistently ignored. Note however, that the axion sector has very interesting 
implications for the strong CP violation problem [l^ . 

^'^As a boundary term, it has no contribution to the equations of motion and will play no role in 
this subsection. However, its contribution is crucial in comparing on-shell actions as we discuss in 
the next subsection. 
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We make the domain-wall ansatz in order to preserve the 4D Lorentz symmetry. 
In the conformal coordinate system, 



ds^ = e2^«(") {dr^ + 7]ijdx'dx^) , $ = $o(r). (3.2) 

Here, r > is the radial coordinate. Boundary is located at r = 0. 

The only non-trivial input in is the dilaton potential V. In order to fix V 
we employ requirements from the dual gauge theory. Holographic dictionary relates 
the scale factor A and the dilaton $ to the energy scale and the 't Hooft coupling 
respectively 

E = e^, Xym = a = e*. (3.3) 
Given these identifications, one can relate the /3-function of the gauge theory to V in 



a one to one fashion[ll|. Although the shape of V{X) is not fixed without knowledge 
of the exact gauge theory /3-function, its UV (small A) and IR (large A) asymptotics 
can be determined. 



UV asymptotics 



In the UV, the input comes from perturbative QCD. We demand asymptotic 
freedom with logarithmic running. This implies in particular that the asymptotic 
UV geometry is that of AdS^ with logarithmic corrections. This requires a (weak- 
coupling) expansion of ^(A) of the form V{X) = 12/£^(H-t>iA + f2A^ + - ■ ■ ). Here i is 
the AdS radius and vt are dimensionless parameters of the potential directly related 
to the perturbative /3-function coefficients of QCD In conformal coordinates, 

close to the AdS^ boundary at r = 0, the metric and dilaton behave as 

where the ellipsis represent higher order corrections that arise from second and 
higher-order terms in the /3-function. The mass scale A is an initial condition for the 
dilaton equation and corresponds to Aqcd just like in the soft-wall model above. 



IR asymptotics 



For any asymptotically AdS space, Einstein's equations dictate the geometry 
in the deep interior be, either another AdS or a singular geometry that terminates 
at r = ro[]T2[. For QCD, it is the second option that is more plausible, as the 
gauge theory is not conformal invariant in the IR. Details of the IR geometry (or 



^^In the latter equation, there is an undetermined proportionahty constant k. However it can be 
set to 1 by a rescaUng in the potential and all physical observables turn out independent of this 
rescaling. Thus, with no loss of generality we can choose Xym = e*. 

^^We will use a "zero" subscript to indicate quantities evaluated at zero temperature. 
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equivalently the large A asymptotics of V) are determined by the requirement of color 
confinement a la [|14|. In particular, we require that the quark-antiquark potential is 
linear. This happens when 



V{X) log^(A), A ^ oo, (3.5) 
and when the parameters Q and P fall into either of the two cases: 



(i) Q > 4/3, P arbitrary — = finite (3.6) 
(ii) Q = 4/3, P > ^ ro = oo. (3.7) 

In the first case, the geometry terminates at finite r, hence this is somewhat similar 
to the hard-wall geometry. The latter case similar to the soft-wall geometry as 
it involves a singularity at r = oo. The asymptotics above also guarantee that 
the classical string configurations do not reach the singularity at vq. The same 
requirement for the particle-like bulk excitations yield an additional condition in 
case (i): Q < 4a/2/3. For Q > 4a/2/3, the glueball spectrum becomes ill-defined[^ 



We note that both the geometry in and in fall into this problematic class. 
These problematic geometries also have undesired features at finite T. 

It is shown in that in both cases, the magnetic quarks are screened and 



the glueball spectrum is gapped and discrete. In case (i) the glueball spectrum 
turns out to be quadratic whereas in case (ii) the spectrum grows as Wn oc n^^. 
For phenomenological reasons, the preferred geometry thus corresponds to the linear 
spectrum with P = 1/2. In this case the asymptotic geometry is, 

dsl ^ e'^'('^)\dr' + dxl) , Ao-e3^/^(^)'Q' (3.8) 



where the constant C is a positive constant related to A in (3.4) 



Parameters of the model 



The dimensionless parameters of the holographic model a priori are (in AdS 
length units): the Planck mass Mpi, which governs the scale of interactions between 
the glueballs in the theory, the scale A£ that plays the role of Aqcd, the string 
length scale is/i and the parameters Vi that specify the shape of the potential V. 
The Planck mass is fixed by studying large T asymptotics, exactly as in eq. ( p.2|) . On 
the other hand, symmetries of the equations guarantee that no physical observable 
depend on Ai. The number i/ig can be determined by comparison with the string 
tension in lattice QCD. For the particular model that is investigated here, this turns 
out to be £/£s ~ 8. This is an encouraging result which shows the a' corrections are 
suppressed by about order 10. 
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Finally, we fix the shape of the potential by arbitrarily picking up a function V 
that satisfies the UV and the IR asymptotics discussed above. A function that does 
the job is, 

V^(A) = ^(^l + viX + V2X^^ log^ (^1 + vsXl + ViX^^^ . (3.9) 

We shall specify the numbers Vi in the following. 

The units in the problem can be fixed by matching the lowest lying 0"''+ glueball 



in our model and in lattice QCD^^ [jT2|. This also fixes the the actual value of the 

2/3 

Planck scale, MpK' . If one wants to compare results of the model with a gauge 
theory with finite N^-, this value gives a cut-off, above which one cannot ignore string 
interactions. For A^^^ = 3 one has MpNc'"^ ^ 2.5 GeV. Of course there is no such a 
cut-off in large QCD. 



4. Thermodynamics of Improved Holographic QCD 

Having defined the theory at zero T, now we look for finite T solutions. At finite 
temperature there exist two distinct types of solutions to the action ( p.l| ) with AdS 
asymptotics, ( |3.4| ): 

i. The thermal graviton gas, obtained by compactifying the Euclidean time in 
the zero temperature solution with r ~ r + 1/T : 

ds^ = e2^o(^) {dr^ + dr^ + dxl) , A = Ao(r). (4.1) 

This solution exists for all T > and corresponds to a confined phase, if the 
gauge theory at zero T confines. 

ii. The black hole (BH) solutions (in Euclidean time) of the form: 

ds' = e^^W (-^ + fir)dT' + dxl^ , A = A(r). (4.2) 

The function /(r) approaches unity close to the boundary at r = 0. There 
exists a singularity in the interior at r = oo that is now hidden by a regular 
horizon at r = where / vanishes. Such solutions correspond to a deconfined 
phase. 

As we discuss below, in confining theories the BH solutions exist only above a 
certain minimum temperature, T > Tmin- 

The thermal gas solution has only two parameters: T and A. The black hole 
solution should also have a similar set of parameters: the equations of motion are 
second order for A and /, and first order for A |^2[. Thus, a priori there are 5 

13 According to ^ this is mo++ 1475 MeV 
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integration constants to be specified. A combination of two integration constants of 
A and A determines A. (The other combination can be removed by reparametrization 
invariance in r). The condition / — >■ 1 on the boundary removes one integration 
constant and demanding regularity at the horizon, r = r^, in the form / — > //i(r/i— r), 
removes another. The remaining integration constant can be taken as fh (or rh, they 
are not independent), related to the temperature by 

4vrT = U. (4.3) 

In the large Nc limit, the saddle point of the action is dominated by one of the 
two types of solutions. In order to determine the one with minimum free energy, we 
need to compare the actions evaluated on solutions i. and ii. with equal temperature. 

We introduce a cutoff boundary at r = e in order to regulate the infinite volume. 



The difference of the two scale factors is given near the boundary as [p^ : 

A{e)-Ao{e) = giT){eA)' + --- (4.4) 
Then the free energy density is given by |T3[: 



Inc = -PNC = lh{M/fK'G{T) - (4.5) 

Here, the entropy density snc is given by the area of the horizon: 

snc = ^n^My^^'^l (4.6) 

One can check (by numerics) that this entropy is precisely the same as follows from 
the 1st law s = —df/dT. This is what one expects as the theory defined by (|3.1j ) 
satisfies the gravitational energy theorems and T is defined in ( |4.3| ) by requiring 
absence of conical singularity at the horizon. 

It is clear from ( |4.5| ) that presence of the first term is crucial for existence of a 
phase transition, as the second term by itself is negative definite. Below, we explain 
the physical meaning of the quantity Q. 



Role of the gluon condensate 

The quantity Q can also be defined from the difference of the dilatons (A = 
exp($)), 

45 

$(6) - $o(e) = -^(r)(eA)Mog(eA) + ■ ■ • (4.7) 

Now, the meaning of Q becomes clear. The AdS / CFT prescription relates bulk fluc- 
tuations with VeVs of dual operators in the gauge theory. As the dilaton couples to 
the operator TrF^, we learn that Q is the difference of VEVs in the gluon condensate 
at finite and zero temperatures: C{T) oc (TrF^)^ — (TrF^)o. 
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Let us perform a consistency check. The dilatation Ward identity in gauge theory 
relates the condensate to the energy-momentum tensor: (T^)to = (TrF^)To- 
The subscript refers to the difference finite and zero T. We shall check this identity 
in the holographic set-up (at leading order in A). The LHS follows from = 
^Nc — '^Pnc- The energy cnc is derived from (p.l|), one finds 

= QO{MpifA^g. (4.8) 

The RHS of the Ward identity is computed by the AdS/CFT prescription: For 
any canonically normalized bulk fluctuation for x(^) = ''^'^^Xoix) + r^+xii^) near 
the boundary, the VeV of the dual operator is {0{x)) = (2A+ — d)xi{x). Taking x 
as 5$ in we find {TtF^)to = 2^iMp^pS^^g. Using the /3-function /3{X) = 

— 6oA^ — . . . we see that this precisely matches the RHS of the Ward identity given 



by 

One may wonder why it works. After all, the prescription is conjectured for 
the pure AdS space and we have a log-corrected AdS here. The reason is that, one 
can generalize the holographic renormalization program of AdS to this geometry by 



explicitly computing the counter-terms |]33[ and show that the contribution from the 



counter-terms cancel out precisely between the finite and zero T components. 

Existence and order of the deconfinement transition 

For a general potential V that obeys the UV and IR asymptotics described in 
the previous subsection, we can prove the following statements: 

i. There exists a phase transition at finite T, if and only if the zero-T theory 



confines as in or ( \3. ID 



ii. This transition is of the first order for all of the confining geometries, with a 
single exception described in Hi: 

iii. In the limit confining geometry P = of l \3. Ao(r) —Cr (as r ^ oo), the 
phase transition is of the second order and happens at T = 3C/Att. 

iv. All of the non-confining geometries at zero T are always in the black hole phase 
at finite T. They exhibit a second order phase transition at T = . 



An heuristic demonstration is given in |T3[ and a general, coordinate independent 



proof will appear in |]32[. Here, let us only mention that the crucial element for the 
phase transition in confining geometries is the existence of (i) a "big" black-hole with 
positive specific heat for small Vh and (ii) a "small" black-hole with negative specific 



One should be careful about the multiplicative factors arising from normalization of $ in (|3 
see il. 
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heat for large r/j. Co-existence of big and small black-holes is just as in AdS BHs 
with spherical horizon. See fig. § for an illustration. It is clear from this figure that 
there exists a for the confining geometries as in eq. ( p.7| , below which both BHs 
disappear. 





Figure 4: Schematic behavior of temperature and free energy as functions of r^, for the 
infinite-r geometries of the type (3.7), for different values of P. 



Numerical Results 



The numerical results that we review in this section are based on ||3^. All the 



thermodynamic properties of the system follow from ( [4.5| ). One numerically solves 
the Einstein-dilaton system for a fixed A^^ and for different r^, corresponding to 
different T (see fig. to obtain snc(T) and G(T). The rest follows from the 
laws of thermodynamics. The potential is chosen in ( |3.9D . Only three of the Vi 
are independent because, as mentioned earlier, the physics is left invariant under 
the rescaling A —>■ kX. We fix one combination of Vi to match the lattice result for 
the latent heat Lh = (0.75Tc)^. The two other parameters are chosen in order to 
obtain good glueball mass ratios^^. A good set of parameters is {f i, f2, "^3, ^4} = 
{0.1, 46, 0.05, 1000}.^^ The rest of the results in this section are predictions. 

We find a transition temperature at ~ 247 MeV which is very close to lattice 
[0. The thermodynamic functions e„c, s„c and pnc are compared with the lattice 
data in fig. ^ left. The temperature dependence of the gluon condensate is shown 
and compared to lattice in fig. ||, right. The speed of sound and the bulk viscosity 
are presented in fig. ^.^^ We conclude that the model presented here is in very good 
agreement with the available lattice data. 



-'^^This is the same in both geometries and fixed by the lowest 0"'""'" mass as A « 290 MeV. 
"'^^We shall not discuss the glueball spectrum here, see |Q. With the potential above, one obtains 
e.g. TOo++*/™o++ = 1-6 which is in well agreement with lattice p5[ |. 

^^The difference in these coefficients and the ones in ^J] are due to a different choice of k here. 
^^This value is for SU{Nc) YM in the large Nc limit which is significantly different from both 



the QCD value and the SW model cited in sec. 2.4. 

^^The derivation of numerical results on the bulk viscosity, along with other dynamical observables 
will appear in [B6|. 
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L/j/r* {Ludni et al. 2005 ) 




5 r. 



Figure 5: Dimensionless thermodynamic functions and the gluon condensate. The dia- 
monds correspond to the lattice data of Boyd et al. 



A last word on the bulk viscosity. Both the low-energy theorems and the lattice 



arguments [|I^ indicate that the bulk viscosity has a peak near Tc. This is what 
we also observe in fig. however the height of the peak is less than the lattice 
evaluation [0^°, see also [|1^ for the same conclusion. 
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Figure 6: Left: Comparison of speed of sound in our model and the lattice result of Boyd 
et al. (diamonds). Right: Comparison of the bulk viscosity wit the lattice data of Meyer. 



5. Discussion and Outlook 

We presented a holographic model for large Nc QCD at finite T, that resolves most of 
the problematic issues of the AdS/QCD models and yields very good agreement with 
the available lattice data. The deconfinement transition results from presence of a 
non-trivial gluon condensate. We also demonstrated that the AdS / CFT prescription 
for computing n-point functions carry over if computed as differences at finite and 

^"Note however that Icomputation of this quantity on the lattice is notoriously difficult and 
afflicted with numerical errors that arise from analytic continuation. 
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zero T. Strictly speaking, the model is valid at large Nc- For finite Nc, there exists 
a UV cut-off, which is about 2.5 GeV for Nc = 3. The a' corrections are somewhat 
under control as the AdS radius is order 10 in string units. However, generally one 
expects corrections from the higher string modes. 

One related problem of all two-derivative effective actions is that the shear vis- 
cosity - entropy ratio is universally fixed as rj/s = 1/Att rather than a function 
of T as expected in QCD. In order to cure this problem, one should consider higher 
derivative corrections in the action. Other possible future directions include study 
of the meson sector via probe D4 branes, turning on a baryon chemical potential 
by charged BHs and eventually searching for explicit non-critical or critical string 
theory backgrounds where the solutions can be embedded. 
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